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The dynamical triangulation model of three-dimensional quantum gravity is shown to have a line of transitions 
' in an expanded phase diagram which includes a coupling /i to the order of the vertices. Monte Carlo renormalization 

, group and finite size scaling techniques are used to locate and characterize this line. Our results indicate that for 

0^ ' n < fii ^ —1.0 the model is always in a crumpled phase independent of the value of the curvature coupling. 

For ^ < the results are in agreement with an approximate mean field treatment. We find evidence that this 
^ . line corresponds to first order transitions extending to positive fi. However the behavior appears to change for 

' fj, > j-i2 2.0 — 4.0. The simplest scenario that is consistent with the data is the existence of a critical end point. 
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I. DYNAMICAL TRIANGULATIONS WITH A MEASURE TERM 



Triangulations provide a discretization of curved Euclidean spacetimes. In the Regge approach, a specific simplicial 
lattice is chosen a priori and the properties of the spacetime are determined by the length of the lattice links. A 
complementary approach, dynamical triangulations, fixes the length of the links to an invariant cut-off and allows 
• the lattice connectivity to determine the geometry of the spacetime. This paper uses the latter approach. A lattice 
representation of the Einstein-Hilbert action is 

S = aNo~ (3Nd (1) 
Oh! 

D , where A^o is the number of vertices and No the number of simplices in the triangulation (the latter is then also the 
i-S^ volume), D is the dimension of the system and a and (i are corresponding chemical potentials. The coupling /3 serves 
*^ \ as the cosmological constant while a is related to Newton's gravitational constant. The connectivity of a triangulation 
plays the same role as the metric in continuous manifolds in the sense that points joined by a link are considered close 
together while those not connected are considered further apart. Continuum theories of quantum gravity generally 
involve an integration over all possible physically incquivalcnt metrics. In dynamical triangulation models, a sum over 
triangulations (i.e a sum over all allowed connectivities) fulfills the same need. Consequently, the partition function 
for the dynamical triangulation model of quantum gravity is 

Z = ^p(r)e^ (2) 

T 

where the sum is over all possible triangulations with fixed topology, S is the action given above and the term p{T) 
allows for the possibility of a non-trivial measure in the space of triangulations. Matter can be coupled to the gravity 
by adding appropriate terms to eqn. (Q) and including a sum over the matter degrees of freedom in eqn. (^) . 

These models of quantum gravity have proven themselves in two dimensions where analytical results are available 
both in the continuum and on the lattice. They agree in all cases where comparisons are possible. This agreement 
occurs with trivial measure p{T) = 1. In higher dimensions, it may be necessary to utilize a non-trivial measure term 
in order to obtain an acceptable lattice theory. In this paper, the three-dimensional model is modified so as to include 
such a measure term. The resulting phase diagram is studied in detail. 

There is additional motivation for adding a measure term. As in lattice gauge theories, it is expected that dynamical 
triangulation models must have a second order phase transition in their phase diagram if they are to have a continuum 
limit and therefore a chance at physical relevance. Many two-dimensional models have the necessary continuous 
transition. However, in three dimensions, while the simple action of eqn. (^ does have a phase transition, it is first 
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order. The phase diagram is one-dimensional (parameterized by a) since f3 is used to fix the volume. Thus, it is 
natural to expand the phase diagram by including some new coupling so that there is a larger territory in which 
to search for continuous phase transitions. Expanded phase diagrams have been examined in the past, produced by 
adding both spin matter Q and gauge matter |^ to the system, but no second order phase transition was found. 
Modifying the action S to include a suitable measure term can be viewed as another attempt to expand the phase 
diagram so as to produce a continuous phase transition. The expectation is that the transition point 04 at /i = will 
extend into a transition line at (/i) in the (a, /i) plane. 

In 1^ a lattice version of a family of measures p{T) is introduced by adding a new term to the action 

S = aNo - PNd + tJ.M (3) 

where 



(4) 



and Oi is the number of D-simplices which include the vertex i, commonly referred to as the "order" of i. This 
additional term in the action corresponds to a continuum measure of the form 



(5) 



being the determinant of the metric). In two dimensions, it is known from Monte Carlo renormalization group 
and exact calculations ||^ that any lattice theory with fi > yields the same continuum limit: this new term 
behaves as an irrelevant operator in the renormalization group sense. In higher dimensions, it is not known whether 
the operator changes the fixed point structure of the theory and if so, how the bare fi should be chosen to approach 
a non-trivial continuum limit. These are the questions we attempt to address here for the case D = 3. 



II. COMPUTATIONAL TECHNIQUES 

We have used Monte Carlo simulations to do computations based on eqn. (3). Three techniques have been most 
useful in analyzing the simulations and drawing conclusions from them: the renormalization group approach, finite 
size scaling, and histograms of the time series data. 

The renormalization group has proven to be a powerful technique for understanding field theory and statistical 
mechanics in flat space. When combined with Monte Carlo simulation it allows a non-perturbative determination of 
the flows and phase structure. A formulation of the renormalization group applicable to dynamical triangulations has 
been developed and established in previous papers It is useful to recall some details of this scheme. One 

iteration of the renormalization group transformation eliminates a single vertex from the lattice. This transformation 
is iterated until the number of vertices in the blocked system is reduced to a target number. Two operators are used 
to monitor the flow in the expectation values: < A^a > and < M >. Their values after blocking depend both on 
the initial volume of the system and the bare couplings. The latter may be used to label different flow lines while 
the former parametrizes distance along such a flow line. For /i < 0, as the renormalization group transformation is 
iterated, all flows initially head toward a common area in the (A'3, Af) plane. Eventually, as the iterations continue, 
the flows diverge and head in directions determined by the phase the couplings correspond to. For large positive /i, 
the flows diverge immediately. 

The renormalization group scheme has been tested previously (in three dimensions) only near the transition a = 
at {^i) for /Lt = 0. It is important to have an independent means of checking whether or not the renormalization 
group prediction of flows, transition couplings, and transition orders is actually correct for /i 7^ 0. If consistent, these 
independent means serve to strengthen the case provided by the renormalization group. Vertex susceptibilities x 
provide just such a check. 

X=[<Nl> - <N^ >^)/N3 (6) 

Close to a phase transition this quantity will typically exhibit a peak. The position of this peak is called the 
pseudocritical coupling. In the infinite volume limit this pseudocritical coupling converges to the true transition 
coupling. The difference between this transition coupling and the pseudocritical coupling is one measure of the 
magnitude of finite size effects. The scaling of the height of the peak in x with system volume can be used to extract 
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a scaling exponent uj. A value for lu of unity would indicate a first order transition. Smaller values are consistent with 
continuous critical behavior. 

If a first order transition is strong enough, it can be seen in the Monte Carlo time series. A double peak structure 
in a histogram of the vertex number data is a classic sign of a first order phase transition. Continuous transitions 
possess only a single peak in this histogram. 

For fixed values of fj,, we search for a transition as a function of a using the above techniques. A variety of values 
of /i are considered, including both positive and negative values. All runs are for (quasi) fixed volumes N3 in the 
range of 4000 to 16000. Typically ten million sweeps were performed where one sweep is defined as A^3 attempted 
elementary moves. For values oi fi ^ A, ten million was found to be much too small a number of sweeps. 

III. RESULTS FOR THE TRANSITION LINE 
A. /i = 0. 

The transition is known to be strongly first order for /i = 0. This is illustrated in fig. Q which shows a Monte 
Carlo time series for the vertex number Nq and in fig. ^ which displays a histogram of the same data. The plots 
show clearly the existence of two metastable states connected by tunneling events typical of a discontinuous phase 
transition. This observation is in agreement with earlier studies |9|. On one side of the transition the system is in a 
crumpled phase and on the other side it is in a branched polymer phase. The crumpled phase is characterized by a 
vertex density < A^q > /A^3 which goes to zero as the volume A3 gets large. In the branched polymer phase this ratio 
approaches 1/3. 
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FIG. 1. Time history of the number of vertices during a simulation at ^ — 0, a = 4.0, and A3 = 8000. Time is in units of 
IOOOA3 attempted updates. 
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FIG. 2. A histogram of the data presented in fig. ^. 

The vertex susceptibility also reveals the first order nature of the transition: in fig.^ results from a system with 
8000 simplices (plotted as circles) clearly display a huge peak. The error bar at the peak is so large because there are 
relatively infrequent tunnelings between metastable states at the transition. Notice that there is a roughly constant 
characteristic value of the susceptibility in the crumpled phase (small a) and a much smaller roughly constant 
characteristic value in the branched polymer phase (large a). These features will be visible at other values of n as 
well. The arrows denote the range of critical coupling determined from a study of the renormalization group flows 
(computed previously |^ and not reproduced here). Notice that it is consistent with the position of the peak in x- 
This is an indication that finite size effects for — 8000 are small at fj, = 0. 

However data from a system with 4000 simplices (with plus as the plotting symbol) is included in the plot for 
comparison. (The error bars are roughly the size of the plotting symbols and were left off for aesthetic purposes). 
The peak is not nearly as high for the smaller system. For significantly smaller systems the peak would become 
indistinguishable from the background. Furthermore, it is clear that the pseudocritical coupling is very far from the 
true critical coupling for this lattice size implying that finite size effects are large. The moral to be drawn from this 
is that it would be very hard to infer the transition order correctly from lattices which are too small. 
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FIG. 3. Vertex susceptibility as a function of a for ^ = 0.0. Arrows indicate the renormalization group prediction for the 
transition region. 



B. /i = -0.50. 

For 7^ the transition value of a has typically not been previously determined (however, see recent work in 
^^ ). The renormalization group flows provide a convenient method of searching for the transition since the flows 
eventual diverge, flowing toward small node density in the crumpled phase and large node density in the branched 
polymer phase. Fig. ^ illustrates this with several values of a chosen to be near the transition. Flows similar to 
the case of /i = are seen here. Each line represents a different value of a and motion along that line (up and to 
the right) represents the evolution of the volume N-j, and the average of M as the initial volume is increased. Since 
the renormalization group is iterated until the number of vertices in the blocked system is the same as the number 
of vertices in the original system when the original volume is 500, systems with larger initial volumes require more 
vertex deletions. Hence, motion along the line corresponds to expectation values for effective theories at larger and 
larger length scales. The values of a are, reading from right to left in the figure, 4.7, 4.8, 4.9, and 5.0. The rightmost 
line, which continues out of the boundaries of the plot, moves toward larger volumes. Since the vertex number is fixed 
this means the vertex density is heading toward zero. This is characteristic of the crumpled phase. The leftmost lines 
cannot go left indefinitely because there is a lower bound on the volume given a number of vertices. These leftmost 
flows correspond to the branched polymer phase. One fiow ends up moving neither right nor left. The coupling for 
this flow is presumably near the transition value. Further points along the flow (i.e. larger volumes) are needed to 
determine which phase it is in or whether it approaches a fixed point. 
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FIG. 4. Renormalization group flows with /i — —0.50. 

At a strongly first order transition, such as in three dimensions with = 0, it is not necessary to go to large 
volumes to see what the phase is, so there are no intermediate flows ||^. In four dimensions, the corresponding 
transition appears second order at small volumes, and there are intermediate flows similar to the one in the figure. At 
larger volumes, there is now evidence that this latter transition is actually first order Thus the three-dimensional 
models at (moderate) negative /x seem to behave somewhat similarly to the the four-dimensional theory at /j, = 0. We 
might expect them to display large finite size effects and this is verified below: it is necessary to go to large volumes 
to get the pseudocritical couplings determined from vertex susceptibilities to agree with the renormalization group 
predictions. 

Fig. H shows the time history for a simulation at the transition and fig. || is a histogram of that data. The existence 
of two metastable states as revealed in these plots lends support to the idea that the transition remains first order as 
fi is varied towards negative values. 
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FIG. 5. Time history of the number of vertices during a simulation at jj. — —0.50, a = 4.85, and — 16000. Time is in 
units of lOOOA^s attempted updates. 
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FIG. 6. A histogram of the data presented in fig. q. 

The vertex susceptibility is shown in fig. ^ There are plateaus of characteristic values in both the crumpled and 
the branched polymer phases. This time, for A'3 — 8000 (with circles as the plotting symbol), the plateaus are 
separated by a much smaller peak. This peak is noticeably to the left of the range of transition values suggested by 
the renormalization group flows (obtained with maximum volume A'3 = 8000). Data for — 16000 (with squares 
as the plotting symbol) show that the transition is moving toward greater a as the volume increases. This is more 
consistent with the renormalization group result. This large change in the pseudocritical coupling is indicative of 
the presence of large finite volume corrections. The large error bar on the peak for A'3 — 16000 is again due to the 
presence of just a few fluctuations between metastable states. 
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FIG. 7. Vertex susceptibility as a function of a for /i — —0.50. Arrows indicate the renormalization group prediction for the 
transition region. 
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C. 



For the limited lattice sizes used in this study the situation becomes somewhat unclear at large negative coupling. 
For moderate lattice volumes there are no clear signals of metastability and a histogram of the vertex number shows 
no strong double peak structure. One might argue that the explanation is that the latent heat is small in this regime 
- perhaps going to zero for some critical /i However, as we have shown our data indicate that finite size effects 
get progressively larger for increasing negative fi. These finite size effects can obscure the true critical behavior - on 
small to moderate lattice sizes the first order transition can appear continuous or even absent. Indeed if we assume 
the correct order parameter for the transition is < Nq/N^ > and if we assume this vanishes in the crumpled phase 
(supported by the mean field arguments we present later) then the observed latent heat is increasing as fi is made 



joint for fi < 0. 

which shows the vertex susceptibility for fi — —0.75. 



more negative. This would seem to argue against a critical end 
Further evidence supporting this scenario is presented in fig. p 
For 8000 simplices (circles), the peak has completely vanished and is replaced with a gradual interpolation between 
the two characteristic values. The range of transition values indicated by the renormalization group is now even 
further to the right of the crossover region in the susceptibilities which is an indication that finite size effects are even 
larger here than for /i — —0.50. The data for 16000 simplices (squares) shows the beginnings of a peak closer to the 
critical regime predicted by the renormalization group method (which used a maximum of 8000 simplices). It appears 
that finite size effects dominate the physics even at this large volume and an even larger volume would be necessary 
to extract unambiguous results. We have not attempted such a simulation. 
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FIG. 8. Vertex susceptibility as a function of a for ^ = —0.75. Arrows indicate the renormalization group prediction for 
the transition region. 



D. = 2.0. 



We have argued that the data favor a first order line for /i < 0. Our evidence concerning the /i dependence of 
finite size effects suggests that simulations for ^ > may suffer smaller finite size effects. This is indeed seen: for 
1.0 < /i < 2.0 it becomes possible to see tunneliiig between metastable states on smaller lattices {N^ = 4000). A time 
history for fi — 2 and — 4000 is given in fig. p| and the corresponding histogram is given in fig. [l^. Similar results 
are obtained at /i = 1 . Analysis of the vertex susceptibility and renormalization group flows also show strong support 
for a first order phase transition at these points on the phase boundary. 
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FIG. 9. Time history of the number of vertices during a simulation at jj. — 2 and = 4000. 
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FIG. 10. A histogram of the data presented in fig. 



E. fj, = 4.0. 



For /Lt = 4, results are quite different. The renormalization group flows are shown in fig. 11. A dashed line has been 
included in the plot as a reminder that the flows all start at N3 = 500. Here the flows diverge immediately rather 
than briefly traveling identical paths as they did for fJ. < 0. The flows fan out more uniformly across the plane. 
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FIG. 11. Renormalization group flows with fi = 4.0. 

Examining the vertex susceptibility (fig. reveals that there is still a peak but that it fails to scale with volume. 
There are no signs of metastability for any value of a. There are also no indications of large finite size effects, so 
the data suggest that there is no transition at this value of /i. Perhaps the first order line terminates for a value of 
fj, between 2 and 4. One other significant feature of the simulations is that there the observed autocorrelation times 
become extremely long. These autocorrelation times could be due to the presence of a nearby critical end point. It is 
also possible that for this value /i we have entered a new phase which is exceptionally difficult to simulate. 
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FIG. 12. Vertex susceptibility for jj. = 4.0 



F. Summary 

The above results as well as some data for other values of /i provide the basis for a plot of the transition through 
the {ct,fi) plane. The results are displayed in fig. The error bars in the figure change direction at a = 9 because 
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for smaller values of a it is convenient to fix \x and vary a when searching for the phase boundary while for larger 
values of a it is more convenient to fix a and vary /i. Notice that the transition appears to head off to large values of 
a as /X is decreased further than /i ^ — 1 . Indeed our data is consistent with a curve which asymptotes to a horizontal 
line. Thus it appears that there is a value of /i below which there is no phase transition as a function of a: there is 
only a crumpled phase. This is because the transition line becomes horizontal. Our results suggest that the transition 
is first order for /i <= 2 with the first order line ending for /i somewhere between 2 and 4. The dotted curve is the 
prediction of our mean field treatment which will be discussed in the following section. 
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FIG. 13. Phase diagram of three-dimensional triangulations with measure term. 

IV. MEAN FIELD ARGUMENTS 

In a previous paper jl^] we have discussed how the two phases of the dynamical triangulation models can be 
understood in terms of a condensation of singular vertices. In general the crumpled phase in D dimensions consists of 
a single (D — 3)-simplex common to a number of Z)-simplices which diverges as a power of the total volume. As the 
vertex coupling a is increased this structure becomes unstable disappearing entirely in the branched polymer phase. 
The transition is driven by the fluctuations in this singular structure |p^ . 

With this as a motivation we have considered an ansatz for the partition function which replaces the counting of 
distinct triangulations with the problem of enumerating the possible orders of the vertices. 



No / No 

rii i—1 \i— 1 / 

The factors uj (rii) give the probability that the vertex i has order rii and the final (5- function constraint reflects the 
fact that the total order is proportional to the total volume |^ Notice that this approach is based on two approximations; 
that the counting of states can be done by counting vertex orders and that the latter can be done in a mean field manner 
by assuming that each varies independently of the others except for the global constraint. Such a model describes well 
certain branched polymer models and has been proposed as relevant to the description of dynamical triangulations 
[15|. This general constrained mean field model was solved in for general probabilities w(ni). As discussed in 
[12|, for describing dynamical triangulations the factors of lo (rii) should be taken to be the number of triangulations 
of the 2-sphere (in three dimensions) surrounding a given vertex i. We take the form [|6l 



Naively c = 4 in D 
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"\^v(i'2{n^~^) (8) 

This has the correct asymptotics {(32 — l^]^^) yields unity for the minimal triangulation of the 2-sphere 
with 4 simplices. Notice that it is now trivial to incorporate the effects of the measure term - the power is simply 
modified from —7/2 to —7/2 + ijl. Using the results of one can write the partition function of this model in the 
thermodynamic limit as 



^ = ^ cxp ((- log A + /32^) cTVa + In (A)iVo) (9) 



where the parameter A is the solution of the equation {p — Nq/N^) 

(10) 



c _ Ai^; (A) 



P (A) 



and, here 



47/2 °° 

q — 2 

This model exhibits two phases; for p < pc = c-pjjj^ the parameter A = 1 and the system is in a collapsed phase 
where vertices of small order behave independently and the global constraint is satisfied by a small number of vertices 
which are shared by a number of simplices on the order of the volume. This regime is identified with the crumpled 
phase of the dynamical triangulation model. Conversely for p > Pc the parameter A varies between zero and one and 
the free energy varies continuously with p. The distribution of vertex orders then behaves as 

-cj(n)e-"'°e^ (12) 

This phase is then identified with the branched polymer phase in the dynamical triangulation model. As p (and 
hence A'o) is varied (by varying the coupling a) the system moves between these two phases. In practice a canonical 
ensemble in which the vertex number can fluctuate is used. 

Zc = Y.e"''"Z{No,N3) (13) 

No 

The phase diagram then contains two phases corresponding to the crumpled and branched polymer phase separated 
by a first order phase transition . In the branched polymer phase A'o /N3 varies continuously with a tending to pc 
as a — at {p) — > 0^. Conversely, Nq/N^ — > in the crumpled phase a < at (p). The phase boundary is predicted to 
be 



at {p)^~ log F,,{1) (14) 

This curve is plotted against the numerically determined phase boundary in fig. |l^ for 1 > p > —1. The agreement 
is rather impressive - at p — the predicted critical coupling is at = 4.1 which is rather close to the best numerical 
determinations at ~ 3.95 — 4.00. It appears for low density (the crumpled phase) the mean field approximation 
captures much of the important physics. The agreement becomes markedly worse for p > I - indeed the mean field 
model would predict that the system was always branched polymer for p > 3/2. This would seem to imply that the 
transition line would once again become horizontal requiring large negative values of a to enter the crumpled phase 
as ^ — > |. We do not see this in the simulations - although the critical behavior of the model certainly seems to 
undergo some form of change for sufficiently large positive p. 

Finally, the mean field phase diagram offers the possibility of understanding the single phase observed at p < pi ^ 
— 1. The critical density pc increases with increasing negative p approaching pc — 1 for p —^ —00. However the vertex 
density cannot increase beyond pmax — ^ - so ii pc > Pmax — 5 the system will always be in a crumpled phase which 
agrees qualitatively with our simulations. 
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V. CONCLUSION 



Our studies of the three-dimensional dynamical triangulation model with measure term allow us to draw several 
conclusions. Firstly, the usual transition point of the model can be extended into a transition line in the {a, /i) 
plane. This line separates a crumpled phase from a branched polymer phase. This boundary corresponds to a line 
of first order transitions. We infer the existence of this line using several methods: renormalization group studies, 
traditional finite size scaling and an approximate mean field treatment. All are consistent. It appears that finite 
volume effects become increasingly important for negative /i. This can make it somewhat problematic to identify a 
first order transition. Indeed the D = 3 simulations for < resemble in this respect the = 4 simulations at 
fi ~ 0. The renormalization group method has shown itself to be extremely useful in this context giving much better 
predictions for such quantities as the infinite volume critical coupling than simple finite size scaling. 

The critical behavior for positive fi appears to undergo a change in the interval /i = 2.0 — 4.0. At ii — 2.0 there 
appears to be a weak first order transition but there is no sign of such a phase transition for fi = 4.0. Nevertheless 
there are very long autocorrelation times there. One simple interpretation of these observations is that the line of 
first order transitions has terminated on a (nearby) critical end point. This might offer the possibility of a continuum 
limit for the model. Other scenarios are possible too. It is possible that the critical line bifurcates into two and a new 
phase appears. Extensive numerical work will be required to distinguish unambiguously between these possibilities. 
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